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ABSTRACT 

We discuss models for the cosmological dark sector in which the energy density of a scalar 
field approximates Einstein's cosmological constant and the scalar field value determines the dark 
matter particle mass by a Yukawa coupling. A model with one dark matter family can be adjusted 
so the observational constraints on the cosmological parameters are close to but different from 
what is predicted by the ACDM model. This may be a useful aid to judging how tightly the 
cosmological parameters are constrained by the new generation of cosmological tests that depend 
on the theory of structure formation. In a model with two families of dark matter particles the 
scalar field may be locked to near zero mass for one family. This can suppress the long-range 
scalar force in the dark sector and eliminate evolution of the effective cosmological constant and 
the mass of the nonrelativistic dark matter particles, making the model close to ACDM, until the 
particle number density becomes low enough to allow the scalar field to evolve. This is a useful 
example of the possibility for complexity in the dark sector. 



Subject headings: cosmology: theory 

1. Introduction 

The striking success of the ACDM model in fit- 
ting the precision WMAP measurements of the 
anisotropy of the 3 K thermal cosmic background 
radiation and the other cosmological tests (Ben- 
nett et al. (2003) and references therein) shows 
this cosmology is a useful approximation to the 
physics of the dark matter and dark energy. How- 
ever, it is not difficult to imagine more complicated 
physics in the dark sector. If the physics differs 
from ACDM enough to matter, it will be manifest 
as anomalies in the fits to the observations. It is 
prudent to anticipate this possibility, by exploring 
models for more complicated physics in the dark 
sector. 

The starting idea for the physics under discus- 
sion in this paper is that the dark matter (DM) 
particle mass may be determined by its interac- 
tion with a scalar field whose energy density is the 
dark energy (DE).^ We explore the physics and as- 

^We do not consider whether such a scenario can emerge 
naturaUy in supersymmetry or axion models. Perhaps it 



trophysics of models for this extension of ACDM 
under the simplifying assumptions of general rel- 
ativity theory, standard physics in the visible sec- 
tor, and, in the dark sector, a Yukawa coupling 
of the DE field to the DM particles. Much of the 
physics, as summarized in §3, is in the literature, 
but we have not seen it all collected and applied 
to the astrophysics. The example application in 
§4, which assumes a single DM family, allows pa- 
rameter choices that make the model predictions 
viable but different from ACDM. The example in 
§5, with two DM families, allows an interesting 
mixture of near equivalence to ACDM at early 
times and complicated departures from this model 
at late times. 

The line of ideas in this topic has a long history, 
which informs assessments of where we are now. 
In §2 we present our selection of the main steps in 
the historical development. 

For convenient reference we write down here the 
forms we will be considering for the action in the 

requires a new view of dark matter. 
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dark sector, as a sum of two terms. The first is 
th.e familiar DE model, 



'DE 



(1) 



where the function V{(l)) of the classical real DE 
scalar field <p is chosen so the field stress-energy 
tensor approximates the effect of the cosmological 
constant A in Einstein's field equation. In numer- 
ical examples we use the power law potential, 



(2) 



where K is a positive constant and the constant 
a may be positive or negative. The DM term, in 
the form used in much of our discussion, is 

Snuf = J d^x^^ [iV^i • d'il) - y{(j) - (f)*)-il>tp] . 

(3) 

The subscript, DMf, indicates this is the action 
written in terms of the wave function ijj for a spin- 
1/2 DM field. In the Yukawa interaction term, y 
is a dimensionless constant and the constant 
has units of energy (with h = 1 = c). If in 
equation (3) is negligibly small, the entire parti- 
cle mass is due to its interaction with the field 
(f). This seems particularly attractive because we 
may need this field anyway, to account for the DE. 
The interaction between the DM and DE allows 
the DM particle mass to be variable, producing 
a long-range nongravitational interaction in the 
dark sector. Both effects can be suppressed by 
the presence of a scond DM family with a differ- 
ent value of (f>^, as we discuss in §5, or by suitable 
choices of 0* and V{(f>) for one family (§4). 

For completeness one might also consider the 
analogue of equation (3) for scalar DM particles. 
When (/)* is negligibly small the analogous form is^ 



'DMb 



= j d'^x^f^g 



(4) 



^The constant in equations (3) and (4) is the sum 
of "bare" and renormalization parts. For generality one 
would add a constant to y^(</) — 4'*)'^ in equation (4). 
Wc ignore considerations of naturalness in the choice of 

these constants. Issues of naturalness and renormalization 
plague other aspects of cosmology and fundamental particle 
physics, as in the meaning of equation (2) within quantum 
field theory, and most notably the value of the vacuum en- 
ergy density. 



In the limiting situation that is thought to 
apply to cosmology, where the DM particle de 
Broglie wavelengths are much smaller than the 
characteristic length scale of variation of the DE 
field, both DM actions (eqs. [3] and [4]) are equiv- 
alent to the model of a classical gas of point-like 
particles with action^ 



4'*\dsi 



(5) 



The invariant interval along the path x^{t) of the 
i**^ particle is dsi = g^,^dxf^dx^ . The DE particle 
mass is rrics = — The absolute value in 
equation (5) has to be a prescription (along with 
y > 0), but as we discuss in §3 it is not needed in 
equations (3) or (4). Equation (5) is a convenient 
form for analyses of structure formation. 

2. Remarks on the History of Ideas 

The particle action in equation (5), with a vari- 
able effective mass, appears in Nordstrom's (1912) 
scalar field model for gravity in Minkowski space- 
time, in the form Lj oc e~'^l'^' dsi (in our notation). 
This form reappears in Misner, Thorne & Wheeler 
(1973), and it still is favored by many, in part be- 
cause the exponential is suggested by superstring 
theory. The functional form does not much affect 
the force generated by the particle-field interac- 
tion, but it affects the cosmic evolution of parti- 
cle masses and the force between them. In our 
preliminary examination of possible alternatives 
to the ACDM model we prefer the linear form in 
equation (5), because it translates to the familiar 
and simple Yukawa interaction in the field action 
model in equation (3). This linear coupling ap- 
pears also (explicitly or as a particular case) in 
many recent discussions of the possible interac- 
tion of DM and DE (e.g. Casas, Garcia-Bellido 
& Quiros 1992; Anderson & Carroll 1997; Bean 
2001, and references therein). 

The particle action with variable mass appears 
in the scalar-tensor gravity theory considered by 
Jordan (1955, 1959) and Brans & Dicke (1961), 
when expressed in units chosen so the action for 
gravity is the Einstein form (Fierz 1956; Dicke 



•'it will be recalled that the exclusion principle affects ini- 
tial conditions - the occupation numbers in single particle 
phase space - but not the equation of motion in the particle 
limit. 
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1965). The units in this theory may be rescaled 
to standard local physics with constant masses 
and a generalized action for gravity, which is the 
form Jordan and Brans & Dicke used to imple- 
ment Dirac's (1938) idea that the strength of the 
gravitational interaction may be small because it 
is rolling toward zero. 

Superstring scenarios led to the thought that 
particle mass ratios (and other dimensionless con- 
stants) may be variable, and in particular grav- 
ity physics and local physics in the visible sector 
may be close to standard while particle masses 
in the dark sector are more significantly variable 
(Damour, Gibbons & Gundlach 1990). This al- 
lows interesting departures from standard cosmol- 
ogy within the tight constraints on gravity physics 
from precision tests in the visible sector. Damour 
et al. (1990) work with a scalar-tensor theory for 
gravity physics, a route taken in many subsequent 
papers. For clarity in our exploratory discussion of 
the dark sector physics we adopt general relativity 
theory, standard physics in the visible sector, and 
a single scalar field that fixes particle masses in 
the dark sector. 

The scalar field model for the DE in equa- 
tion (1) was introduced by Wetterich (1988) and 
Peebles & Ratra (1988) (as reviewed in Peebles & 
Ratra 2003). This DE model allows one to imagine 
that the effective cosmological constant is evolv- 
ing to its "natural" value, A = 0, and is small 
now because the universe is old, a natural exten- 
sion of the ideas of Dirac (1938), Jordan (1955), 
and Dicke (1964). Wetterich (1995) seems to have 
been the first to propose that the scalar field in 
this model for the DE may also fix the DM par- 
ticle mass. The idea has since been discussed in 
a considerable variety of contexts (e.g. Damour, 
Piazza & Vencziano 2002, references therein, and 
the references in the following discussion). 

A scalar interaction present only in the dark 
sector makes the accelerations of visible and dark 
matter test particles different, an effect we will 
call a "fifth force." The empirical constraints on 
this kind of fifth force arc considerably weaker 
than the constraints from the Eotvos experiment 
in the visible sector, as was recognized from the 
beginning of the modern discussions, in Damour 
et al. (1990). The first numerical example we 
have seen of the effect of this kind of fifth force 
on the growth of mass density fluctuations in the 



expanding universe is in Amendola (2000). Amcn- 
dola & Tocchini-Valentini (2002) point out that 
the fifth force in the dark sector might have a 
substantial effect on the relative distributions of 
baryonic and dark matter. But we now have con- 
vincing evidence (Bennett et al. 2003 and refer- 
ences therein) that structure grew out of primeval 
adiabatic departures from homogeneity, and good 
evidence that the growth of the mass density fluc- 
tuations is not very different from the ACDM pre- 
diction, from the consistency within this model 
between the power spectra of the present distri- 
butions of galaxies and the 3 K thermal cosmic 
background radiation (the GBR). Thus one is in- 
terested in DM-DE interaction models that can be 
adjusted so the fifth force and the evolution of the 
DE field are weak enough to fit the now demand- 
ing observational constraints, but strong enough 
to make an observationally interesting departure 
from AGDM. 

The search for models has been influenced by 
the attractor concept, that the physics may have 
the property that the astrophysics is insensitive to 
initial conditions. Peebles & Ratra (1988) intro- 
duced the DE power law potential in equation (2) 
because it has this attractor property. The physics 
need not have an attractor, of course, as in the 
example of Franga & Rosenfeld (2002), which is 
based on Wettrich's (1988) potential and suitably 
chosen initial conditions. We discuss models of 
this kind in §4. 

In the attractor model considered by Anderson 
& GarroU (1997), and more recently by Comelli, 
Pietroni & Riotto (2003), the potential of the DE 
field is the sum of a term linear in the field and 
proportional to the DM particle number density 
(as in equations [3] and [5]) and a power law self- 
interaction term with a > in equation (2). The 
field is assumed to have been attracted to the min- 
imum of the total potential. Within the models 
considered in §4 this attractor case is unaccept- 
able, because the fifth force is too large. 

Attractor solutions may also be relevant to the 
dilaton potential. In the scenario considered by 
Damour & Polyakov (1994) the masses of all par- 
ticles have minima as a function of the dilaton 
field at a universal value near this minimum 
the fifth force scales as (</> — (/'m)^- They show that 
comic evolution can cause the present value (j) of 
the dilaton field to be close to thus suppress- 
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ing the fifth force. 

In the class of models we are considering, the 
contribution to the potential energy of the DE 
coming from its interaction with the DM has a 
minimum at zero DM particle mass. This mini- 
mum value certainly is not acceptable for the DM 
in galaxies, but one can imagine that there are 
two families of DM particles, with different values 
of (/)* in equations (3) and (5). At largo enough 
DM particle number densities the DE field would 
be locked to the zero of the particle mass for the 
preponderant family, making these particles rela- 
tivistic (and with a mass density that can be ac- 
ceptably small). This has the effect of suppressing 
the fifth force and eliminating the evolution of the 
massive DM particle mass and the evolution of the 
DE density. This resembles Damour & Polyakov's 
(1994) "least coupling," but with the difference 
that in an expanding universe the particle num- 
ber density must eventually become low enough 
to release the field. 

We see in this history of ideas a conservative 
aspect of theoretical physics. The Nordstrom ac- 
tion, embodying a variable particle mass, was un- 
der discussion before Einstein had completed his 
general relativity theory of gravity. It reappeared 
in the 1950s and 1960s, in scalar-tensor general- 
izations of general relativity expressed in terms 
of the Einstein action. The scalar-tensor theories 
were developed to explore Dirac's idea, that the 
strength of the gravitational interaction may be 
variable, and these theories served also as a guide 
to the work of developing precision tests of gravity 
physics. The current reappearance of this action 
is motivated in part by superstring scenarios and 
in part by the continuing fascination with variable 
parameters of Nature. Moreover, the Nordstrom 
action in the form of equation (5) is equivalent to 
a Yukawa interaction with a classical scalar field, 
the form of which was introduced for very different 
purposes in meson and weak interaction physics. 
To be discovered is whether whether this conver- 
gence of ideas from particle and gravity physics is a 
useful guide to observationally significant aspects 
of the physics of the dark sector. 

3. Basic Relations 

We begin with the physics of the DM particle 
model in equation (5). Many of the results sum- 



marized here have appeared in one or more of the 
papers cited above, but we have not seen them all 
collected or applied. The relation of the particle 
model to the field model in equation (3), which is 
discussed in §3.4, might be considered self-evident, 
but it should be checked in the present context. 

3.1. Particle and field equations 

We simplify notation in this subsection by set- 
ting to zero (which has the effect of shifting 
the minimum of ¥{(())) and taking <j) to be posi- 
tive. Throughout y is positive. 

The particle action in equation (5) gives the 
equation of motion 

d dx" y(p dgp„ dx^ dx'^ dycp 

^#ff^.^ = Ya^^^ + a^- (6) 

We leave the constant y in this equation because it 
is useful to note that the particle four-momentum 
is = y^adx^^/ds. When spacetime curvature 
fluctuations can be neglected the equation of mo- 
tion is 

dap d_ ayct^^ yiT^g^ (7) 

where a{t) is the cosmological expansion factor as 
a function of the proper world time t, and the 
proper peculiar velocity is v = adx/dt. When 
the spatial variation of the DE field cj) may be ne- 
glected the momentum is conserved; if in addition 
the proper peculiar velocity is nonrelativistic, the 
velocity scales as u cx l/{a{t)4>{t)). 

The dark matter that is bound to galaxies and 
clusters of galaxies has to be nonrelativistic. These 
systems are well described by the weak field limit 
of gravity, where the gravitational potential satis- 
fies 

V2<I./a2 = ATrGpf,{t)S{x, t). (8) 

The mean (background) nonrelativistic mass den- 
sity is Pb{t), and S = Sp/ph is the mass density 
contrast. We shall write the DE field as 

(t> = Mt) + M^,t), (9) 

where the mean background field 4>b{t) is a func- 
tion of world time t, and the departure 01 from ho- 
mogeneity may be treated in linear perturbation 
theory. (We will also be assuming \ <C \^b — ^*\, 
so that there is no risk that the field value passes 
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through zero, which would make the DM tran- 
siently relativistic.) In linear theory we may also 
neglect the term proportional to v-Vcj) in equation 
(6). With all these approximations equation (6) 

becomes 



dv 
It 



a (t>b 



a V 4>b{t) 



. (10) 



where the dot is the derivative with respect to 
time t. The expansion of the universe produces 
the familiar slowing of the peculiar velocity in the 
second term of this equation, while the evolution 
of the DE field value produces a term that may 
increase or decrease the peculiar velocities. The 
spatial variation of the DE field produces the fifth 
force term in the right hand side of the equation. 
This force tends to move DM particles so as to 
minimize their masses y(j). 

The DE field equation from the action in equa- 
tions (1) and (5) is 

1 



d , d(\> dV dVi „ 



The term from the interaction with the DM is 



dVi 

^4 



= vE 



dsi (5(x — Xi) 
dt ,7=0 



y 



J2Jl-vM^-^0/a'. (12) 



The last expression neglects the effect of spacetime 
curvature fluctuations on the DM source term for 
the DE. 

One sees that when the particles are relativistic, 
Vi ^ 1, the source term dVi / d(j} vanishes. Damour 
& Polyakov (1994) express this in terms of an 
equation of state. We find it convenient to use 
instead the proper DM particle number density, 

n(x, t)=Y, '5(x - Xi)/a3 = ^ 5(r - n), (13) 

i i 

where 6y — a(5x is a proper relative position. Thus 
the source term may be written 



Z = yn{yi,t){\/l- 



(14) 



The brackets signify the mean of the reciprocal 
Lorentz factor, 7, for the DM particles. The 



pressure in this homogeneous gas of DM parti- 
cles is p = ^y(pnv'^/3 and the energy density is 
p = jy(f>n, so another form for the source term is 
dVi/dcj) — {p— 'ip)/4i- This is the form derived by 
Damour & Polyakov (1994). 

3.2. The Fifth Force 

There is a fifth force in the dark sector, which 
we analyze for a single nonrelativistic DM family 
model. The DM in galaxies and clusters of galaxies 
is well described by the approximations of equa- 
tions (8) and (10), which ignore relativistic corrc- 
tions and the effect of spacetime curvature fluctu- 
ations on the DE field equation. Returning now 
to the general case when is non-zero, we can 
write the spatial mean of the field equation as 

where (j)b is the mean field (eq. [9]) and nb{t) is the 
mean particle number density. Here and below, 
the last term has the sign of — 0* . The departure 
(1)1 from the homogeneous part of the fleld satisfles 
the equation 



df^ 



1 



a dt 



d^V 



±ynbS = 0. 

(16) 

The DM number density contrast 5n/n has been 
replaced by the mass contrast 6 = 5p/ p, because, 
as we argue below, in situations of interest the 
fractional perturbation to the fleld (j) is small com- 
pared to 5n/n. 

For the analysis of structure formation at mod- 
est redshifts we are interested in density fluctu- 
ations on scales small compared to the Hubble 
length, which means the time derivatives in equa- 
tion (16) are small compared to the space deriva- 
tives. The term d'^V/d(jP' is also small in many 
cases of interest, so a useful approximation to the 
field equation is 



(17) 



It follows from equations (8), (10) and (17) with 
Pb = y\(l>b — 4>*Wb that the ratio of the fifth force 
to the gravitational force in the dark sector is 



(3: 



|V01 



|(/)6-<^.|V$ 47rG((/.6-<^*)2- 



(18) 



5 



The evolution of the mass density contrast in lin- 
ear perturbation theory satisfies 



85 



a {(j)b - 0,) / dt 



= 4TrGpb{l + P)6. 

(19) 

This differs from the usual expression by the fac- 
tor 1 + (3, that takes account of the fifth force, 
and the part cpi, / {(p^ — 0, ) , from the evolving DM 
particle mass (as has been discussed by Amendola 
& Tocchini-Valentini 2002; Matarrese, Pietroni & 
Schimd 2003, and others). 

Now we can check that 5n/n ~ 5p/ p. Equa- 
tion (17) applied to a particle concentration with 
contrast 5 and size r says 



5m (pi —ynbr'^5 

m {4>b-4>*) \4'b-<t>*\' 



(20) 



where H is the Hubble parameter. The sign of 
the mass shift 5m /m is opposite to the sign of 5, 
consistent with the attractive nature of the fifth 
force. Using the condition that the DM mass den- 
sity is not greater than the total, we find that the 
fractional mass shift is 



< 



l3{Hrf\5\- 



(21) 



This is small because (3 cannot be much larger 
than unity, and we are interested in density fluc- 
tuations on scales small compared to the Hubble 
length. 

3.3. A Relativistic Dark Matter Family 

The DM bound to galaxies and clusters of 
galaxies has to be nonrelativistic, with (Vl — v"^) 
close to unity, so the source term (14) due to the 
scalar field interaction with such DM is simply 
proportional to the particle number density. But 
there may be more than one DM family, and (j) may 
be drawn close to the zero of the source term be- 
longing to one of the families, causing that family 
to be relativistic. Here we consider the behavior 
of this piece of dVi / d(p function of the field 
value, under the simplifying assumption that the 
scalar field is spatially homogeneous. In §5, on a 
model with two DM families, we consider the re- 
sponse of (j) to the source term and the effects of 
inhomogeneities . 



The velocity of a DM particle in this relativistic 
family is 

k/a 

V = ' =, (22) 

where the cornoving wavenumber of a DM parti- 
cle is k and the proper peculiar momentum is kja. 
The initial distribution of comoving wavenumbers 
is determined by the DM particle production pro- 
cess. We will discuss the evolution of ^ in §5. For 
the purpose of this subsection we are concerned 
with the case y\(p\ < k/a <C H^^. The first in- 
equality means the DM particles are relativistic; 
the second means the discussion of the evolution 
of the mass distribution can ignore time deriva- 
tives of (j) and the expansion of the universe. In 
this approximation the energy of a DM particle 
is conserved as it propagates through space and 
encounters DE field gradients. Therefore the lo- 
cal value of the DM particle velocity changes in 
response to changes in the local value of the DM 
particle mass, as determined by the local value of 
the DE field. With this in mind, we can rewrite 
equation (14) as 

^=yn{t){yHy'^' + kyaT'^^),, (23) 

where the average is over the distribution of parti- 
cle wavenumbers in the relativistic family. When 
(j) passes through zero (or more generally, through 
the zero of the particle mass) in a close to homo- 
geneous way it does not greatly aflFect the distri- 
bution of momenta k/a, but it makes the DM pe- 
culiar motions relativistic. That causes the source 
term (23) to vary smoothly from dVj /d(j) = yn at 
ycj) » kes/a to dVi/d(p = —yn at ycj) <C —kes/a, 
where the effective comoving wavenumber fcefr is 
defined by the average in equation (23). 

As a final remark, we note that, according to 
the action principle, the equation of motion is de- 
rived by extremizing the action (4) with the par- 
ticle orbits held fixed. Thus although Vj may be 
written as a function of (j), n, and the momenta 
k/a, the source term dVi/d(j) is not the same as 
the derivative of Vj with respect to <p at fixed n 
and k/a. In the derivation of the field equations 
from the Lagrangian (3), which will be discussed 
in §3.4, the independent field V is held fixed as (p 
is varied, leading to equation (14). 
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3.4. The field action 

Here wc consider, in the limit where the parti- 
cle de Broglie wavelengths are small compared to 
the length scale of variation of (j), the equivalence 
of the particle action in equation (5) to the spin- 
1/2 fermion field action in equation (3) and the 
boson field action in equation (4), for the purpose 
of deriving particle orbits and the source term for 
(j). For brevity we again take (j) to be positive and 
suppress 0*. 

If the parameters y and </> in the Yukawa inter- 
action term in equation (3) are positive the parti- 
cles created by the field ijj manifestly have positive 
mass y(j). A chiral rotation by tt changes the sign of 
V'')/' without changing the kinetic part of the field 
equation, so when ycj) is negative the chiral rota- 
tion yields the usual sign for the mass in the Dirac 
equation. Alternatively, we can leave the "wrong" 
sign for the mass when 4> is negative, and note by 
the chiral transformation argument that the solu- 
tions to the field equation in this case make V'V' 
negative, meaning y(j)'4)'tp is never negative. This 
condition leads to the prescription for the absolute 
field value in the particle action (5). 

In the field action (eq. [3]) the source term for 
the DE field is 

dVi/d(j) = ytptp = y n\/l - v^, (24) 

when y and (j) are positive. The last expression 
follows because ^jxJj and — both arc scalars, 
and we know they are the same - up to the sign - 
in the nonrelativistic limit. ^ This agrees with the 
source term in equation (14) in the particle model. 

The second step in demonstrating the equiva- 
lence of the actions (3) and (5) is to check the 
equation of motion of wave packets. An easy way 
to proceed uses the commutator of the particle 
momentum operator with the Dirac Hamiltonian 
H = a-p + Pycj), 

[p,H] = -iypW<p. (25) 

It follows that the time derivative of the expecta- 
tion value of the momentum is 

^ = -yVcj) I Sr^^ji^l). (26) 
at J 

^Onc can check these arguments by writing down plane wave 
solutions to the Dirac equation. The DE source term in 
equation (24) can be derived from the free quantum field 
operator for i/i, apart from the standard problem with the 
zero-point contribution to the particle number operator. 



But the last factor is the integral over ijjijj, which 
we know is the reciprocal of the Lorentz fac- 
tor for a single particle wave packet, as in equa- 
tion (24). Equation (26) thus agrees with the rate 
of change of momentum in the particle model in 
equation (7). 

For completeness let us check the relation be- 
tween the momentum and velocity of the wave 
packet. An easy way is to use a WKB approx- 
imation. In considering the motion of the wave 
packet we can ignore the time evolution of (j). The 
interesting part of the spatial variation in one di- 
mension of a wave function with energy e is 

^-^-^^{f (27) 

which means the momentum defined by the gradi- 
ent operator is 

p = - (28) 

as usual. The velocity of a wave packet con- 
structed as a linear combination of these energy 
eigenstates follows from the stationary point of 
the exponential: v = p/e, again as usual. And 
these two results give the standard relation be- 
tween momentum and velocity of a particle of mass 
ycj). On can also check that, when the time evo- 
lution of (j) and the expansion parameter a{t) can 
be neglected, the time derivative of equation (28) 
agrees with the rate of change of momentum in 
equation (7). 

Similar arguments show that in the limit of 

short de Broglie wavelengths the boson DM model 
in equation (4) also reduces to the point particle 
model (apart from the problem of the zero point 
contributions to 0^ and y^), with DM particle 

mass y\(j) — 4>*\- 

4. Models with one dark matter family 

In the models presented in this section, the DE 
potential is the power law form in equation (2), 
with positive or negative power law index a, and 
there is one family of nonrelativistic DM particles 
with = 0, meaning the DM particle mass is 
y times the DE field value. For the purpose of 
this preliminary exploration we neglect the mass 
in baryons, so the matter density parameter is 

^mHl = ^TTGpbito) = ^nGyMtohbito). (29) 
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Here and below the subscript o means the present 
value. The Hubble parameter at the present world 
time to is Ho- Throughout we use flm = 0.3 and 
we assume space curvature vanishes. ^ 

We use the dimensionless variables 



r = Hot, f = G^/^^b, 



(30) 



in terms of which equation (15) for the mean field 
h is 

(Pf 3 da df aK 3 /ao\3 17, 



dr"^ adrdT Stt 



(^) (31) 
V o / fo 



and the Friedman equation for the expansion rate 
is 



adr J 



fo^aJ ^3 



2\Jt) +7^ 

(32) 

where the dimensionless parameter representing 
the constant K in the power law potential V{<j)) is 



(33) 



The present values fo and Go of the field and 
expansion parameter appear in the combination 
fo/al, which is the unknown final condition for 
given initial conditions. 

At high redshift the first term on the right hand 
side of equation (31), which represents dV/dcf), is 
relatively small. When this term may be neglected 
the first integral of equation (31) is, apart from the 
decaying term, 



dT 



Stt/, 



'1 (^Y 

o ^ 0"' 



(34) 



At 2; > Zeq, where the expansion is dominated by 
radiation, the expansion factor varies a oc r^/^, 
and equation (34) says the departure from the 
initial value of / grows as r^/^. At lower red- 
shift where the expansion is matter- dominated, 



®For simplicity in this exploratory discussion of observa- 
tional constraints we do not adjust the value of Qm to take 
account of the fifth force. In a theory with a fifth force 
/? enters different measures of Qm in different ways. For 
instance, the relative motions of dark matter halos depend 
on the product Clm{l whereas weak lensing depends 

on the fifth force only indirectly, through whatever efl^ect 
the scalar field has on the angular size distance. The dy- 
namics of ordinary matter are not directly affected by the 
fifth force in the dark sector. 



a oc r^/'"', the departure grows as logr in the ap- 
proximation of equation (34). At still lower red- 
shifts dV/d^ may be important, and we need a 
numerical solution. 

We commence the numerical solution at a fixed 
initial time corresponding to equality of mass den- 
sities in matter and radiation in the ACDM model. 
We start with an arbitrary choice for the initial 
field value /, = G^/^^,. The initial value of df /dr 
is taken from equation (34). The final field value 
(and hence energy density) has to be consistent 
with da/dr = a at the present epoch. We achieve 
this by iteratively adjusting k. (We found this 
more convenient than choosing k. and seeking the 
initial field value.) Having solved numerically for 
/(r), we can find the epoch Zeq at equal mass den- 
sities in matter and radiation. 

Table 1 lists parameters and present values of 
some quantities of interest for solutions with three 
choices for the value of a, omitting numbers that 
are so far off the ACDM model prediction as to 
seem uninteresting. The second column is the ini- 
tial field value, expressed in units of the Planck 
mass. The third column is the value of k required 
for a consistent solution. The fourth column is the 
ratio of the field value at Zcq to the present value. 
Because we are assuming = 0, this is the ratio 
of DM particle masses then and now. The redshift 
at equality scales in proportion to this ratio. 

Figures 1 and 2 show the evolution of the DE 
field in solutions with a = — 2 and a = 6. The lat- 
ter look much like the solutions for a = A entered 
in Table 1. Since we have set = 0, the DE field 
in solutions with a < is drawn toward zero. At 
a = —2 and the smallest initial field value listed 
in the first line of the table and shown as the dot- 
ted curve in Fig. 1 - the field has passed through 
zero slightly before the present epoch. Among 
other undesirable consequences, this would have 
made the DM transiently relativistic, driving the 
DM out of the halos of galaxies. The slightly larger 
initial field value in the second line of the table, 
with the appropriate adjustment of k, removes this 
problem: in this case, shown as the dashed curve 
in Figure 1, the DM particle mass has not yet 
passed through zero, but that will soon happen 
and the halos will be disrupted. 

When a is positive the potential for cp has a 
minimum away from zero. If the initial value of 
(p is small enough the field relaxes to this mini- 
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Table 1 

Numerical results for one DM family. 



a. 




K 




X 


'peak 




<5acdm 


(acdm 


-2 


0.95 


1.5E+01 








-2 


1.00 


1.4E+00 








-2 


2.00 


3.0E-02 


1.22 


1.18 


1.18 


-2 


4.00 


5.7E-03 


1.04 


1.04 


1.04 


4 


0.50 


4.7E-03 








4 


1.00 


3.5E-03 


2.17 


2.88 


2.02 


4 


2.00 


6.8E-01 


1.19 


1.18 


1.17 


4 


4.00 


1.8E+01 


1.04 


1.04 


1.04 


6 


0.50 


l.OE-02 








6 


1.00 


1.3E-02 


1.33 


1.37 


1.29 


6 


2.00 


2.0E+00 


1.18 


1.17 


1.16 


6 


4.00 


2.8E+02 


1.04 


1.03 


1.04 



mum by the present epoch. This is seen in our 
solution with the smallest initial field value, listed 
in the fourth entry from the bottom of Table 1 
and plotted as the dotted curve in Fig. 2. In 
this case the field oscillates about and approaches 
the minimum of the potential. This solution is un- 
acceptable, however, because the relatively small 
field values produce a large fifth force on the dark 
matter, which substantially enhances the growth 
of mass density fluctuations, as we discuss next. 
At the two largest initial field values in Table 1, 
the solutions for a = 6 are well away from the 
minimum of the potential. They look much like 
the solutions for a = —2, and as we show next, 
produce only a modest effect on the evolution of 
mass density fluctuations.^ 

Figures 3 and 4 show numerical solutions to 

'Anderson & Carroll (1997) and Comelli, Pietroni & Riotto 
(2003) consider the case that the scalar field sits at the 
minimum of the effective potential. They adopt the same 
particle coupling and potential (with a > 0) as in the ex- 
amples in this section. However when the scalar field is at 
the minimum of the effective potential — as in an attractor 
scenario - there is an unacceptably rapid evolution of the 
DM particle mass subsequent to decoupling. This evolu- 
tion, and the fifth force, can be suppressed by choosing a 
large value of </>« ~ mpj in equation (18). Our acceptable- 
looking cases arc not attractor solutions: tlicy arc sensitive 
to the initial value of 0, and their success relies in part on 
the assumption that (j> is far from its value at the minimum 
of the potential. 



equation (19) for the evolution of the mass density 
contrast 5{t) in linear perturbation theory. The so- 
lutions arc normalized to a common initial value 
at rcdshift z = 1300, roughly the epoch of decou- 
pling, and they are multiplied by the redshift fac- 
tor 1 + z = ao/a{t) to scale out the main trend of 
the evolution. The solution for the ACDM model, 
with the same value of flm, is plotted as the short 
dashed curves in both flgures. The solution for 
a = -2 and ft = G'^/'^(f>i = 0.95 is not plotted 
because equation (19) does not take account of 
the transient rclativistic motions of the DM par- 
ticles. The flfth column of Table 1 lists the ratio 
of the growth factor since decoupling, 6o/Sdec, to 
the prediction of the ACDM model. 

At the two largest initial flcld values and all 
three choices of a, the growth of density fluctu- 
ations is close to the ACDM prediction. In the 
solution for fi = I and a = —2, plotted as the 
long dashed curves in Figs. 1 and 3, the growth 
of the density contrast is about four times that 
of the ACDM model. That is ruled out by the 
consistency of the CBR temperature anisotropy 
and the large-scale fluctuations in the galaxy dis- 
tribution within the ACDM model. For = 1 
and a = 6 the density fluctuation growth factor 
is more than a factor of two different from ACDM 
at redshift z = 10, but happens to be fairly close 
at the present epoch (as one sees by comparing 
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the long and short dashed curves in Fig. 4). This 
solution is challenged by the position of the peak 
of the CBR fluctuation spectrum, however, as we 
now discuss. 

We estimate the angular scale of the peak of 
the CBR temperature anisotropy power spectrum 
as follows. Since the physical wavelength of the 
mode that produces the peak of the fluctuation 
spectrum is set by the Hubble length at Zeq, the 
wavenumber at the peak varies with the model 
parameters as 



k 



peak 



ao t 



eq 



do f, 

OC CX — 

^cq Jo 



(35) 



The second step follows because the expansion 
time is inversely proportional to the square of the 
temperature at Zcq, that is, icq oc a^q, and the last 
step follows because the redshift at equal mass 
densities in radiation and DM varies as /eq//o 
through the evolution of the DM particle mass. 
The peak of the angular power spectrum of the 
CBR temperature is at spherical harmonic num- 
ber Ipcak ~ kpcakf, whcrc the angular size distance 
r = J dt/ a is integrated from decoupling to the 
present epoch. Thus the ratio of the spherical har- 
monic index Zpcak in the model to the predicted 
index at the peak in the ACDM model with the 
same cosmological parameters is 



peak 



fet 



^ACDM fo J'ACDM 



(36) 



This ratio is listed in the last column of Table 1 . 

The model with a = 6 and G^^^^i = 1, whose 
present density fluctuations happen to be close to 
the ACDM prediction, puts the peak of the CBR 
temperature fluctuation spectrum at angular scale 
~ 30% smaller than ACDM, which likely is unac- 
ceptable. At G^^'^(j)i — 2 the peak is shifted from 
ACDM by about 16% , which may be tolerable 
within the uncertainties allowed by the other cos- 
mological parameters that determine the value of 
^peak- A closer analysis of the joint distribution 
of allowed values of (pi and the cosmological pa- 
rameters seems inappropriate in this preliminary 
exploration. The point we wish to demonstrate is 
that there is a range of initial field values that pro- 
duce a significant but acceptable departure from 
the behavior of the ACDM model. 

The evolution of the DE density nowadays is 
characterized by an effective DE equation of state. 
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0.001 0.01 0.1 

expansion parameter a(t) 



Fig. 1. — Evolution of the DE field in models with 
the power law exponent a = — 2 in the potential 
(eq. [2]). The solutions are fixed by the initial 
value of / = G^^'^cj) listed in the first four rows of 
Table 1. The initial values are close to the field 
values at the left-hand edge of the plot. 




0.001 



0.01 0.1 
expansion parameter a(t) 



Fig. 2. — The same as Fig. 1 for solutions with 
a = 6. The initial values of (p are listed in the last 
four entries in Table 1. 
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expansion parameter a(t) 



Fig. 3. — Evolution of the mass density contrast in 
linear perturbation theory in models with a = —2. 
The density contrast has been multiplied by the 
redshift factor 1 + z to scale out the evolution 
when the expansion is matter-dominated. The 
short dashed curve is the solution for the ACDM 
model. The line types of the other curves match 
Fig. 1, where the initial field values are close to 
what is plotted at the left side of the figure. 



' \ 




0.001 0.01 0.1 

expansion parameter a(t) 



Fig. 4. — The same as Fig. 3 for solutions with 
a = 6 



Fig. 5. — The evolution of the equation of state 
parameter (Fig. [38]) in solutions with a = —2. 
The line types match Fig. 1. 



Q. 




expansion parameter a(t) 



Fig. 6. — The evolution of the equation of state 
parameter in solutions with a = 6. The line types 
match Fig. 2. 
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We define the effective pressure pcs by the expres- 
sion for local energy conservation, 

^(Pdm + Pde) = -3-(pdm + Pde +Pefr)- (37) 

The ratio of the effective pressure to the DE den- 
sity is 

^^p^^_y-0V2+(V0)V2_ (38) 
Pde F-F0V2-F (V0)2/2 

Figs. 5 and 6 show the evolution of the equation 

of state parameter w in our numerical solutions, 
which neglect the gradient energy density in equa- 
tion (38). The complicated behavior of w in the 
solutions with the two smallest initial field values 
is of no interest because the models are not viable. 
At the two larger initial field values the parame- 
ter is close to constant at w ~ — 1 in the range of 
redshifts reached by the SNela observations. At 
higher redshifts w ~ because in these solu- 
tions the DE energy is dominated by (p^ /2, but at 
high redshift the DE density is well below the DM 
mass density. 

5. Cosmologies with two dark matter fam- 
ilies 

Here our field Lagrangian is 

+i-ipj ■ dip + iipsj ■ dips 

-y{(p*- 4>)i>ip -ys(p'4)sips-, (39) 

where for definiteness we choose y,ys,(ps and <?!)« to 
be positive constants. The mean number densities 
are n in the first family, with wave functions ip, 
and Us in the second family. If (p is never close to 
zero or to <p<, then there are two nonrelativistic DM 
families, a situation that can be equivalent to what 
we discussed in the previous section. A new possi- 
bility offered by this model is that at high particle 
number density the DE field is "locked" to (p — ^ 
to minimize the effective potential of the family 
with the larger value of the Yukawa coupling con- 
stant times the number density, making that fam- 
ily relativistic. If ygUs > yn in equation (39), and 
V{(p) is subdominant, the particle mass in the mas- 
sive DM family is fixed to m = y(p^, and the DE 
density is fixed to puE = Vi^Ps), so the DE behaves 
like Einstein's cosmological constant. We show in 



§5.1 that the locking can also substantially sup- 
press the fifth force. In §5.2 we comment on the 
complicated behavior when the number density in 
the second family becomes small enough to allow 
the DE field to evolve. 

5.1. A relativistic family in the dark sector 

We make several approximations that simplify 
the analysis of the fifth force when the DE field 
is close to the minimum of the potential of one 
family. 

First, we assume the length scale of the density 
fluctuations of interest is much smaller than the 
Hubble length, so the dark matter particles can 
relax to near time-independent equilibrium. Sec- 
ond, we assume that at high redshift, when the 
DM number densities arc large, the scalar field re- 
laxes to the minimum of the particle potential. We 
take ygUs '> yn, so the mean value of the scalar 
field is close to zero. This means there is a nonrel- 
ativistic DM family with mean mass density 2/0* n 
together with a relativistic second family. Third, 
it is reasonable to take it that the space distri- 
bution of the relativistic family is close to homo- 
geneous. We can see how this comes about by 
considering the distribution of positions and mo- 
menta p = — 2/2(/>(x)2 in single particle phase 
space. If the distribution has relaxed to become 
nearly independent of time and a function only of 
the energy e, then the space number density distri- 
bution for the second family particles with energy 
e is 

ns(x) (X j d^p6{e - + p^) x p/e = v. (40) 

This must be averaged over the distribution of en- 
ergies. But we see the key and familiar point, that 
when the second family is relativistic, that is, v is 
close to unity, the space distribution is close to 
homogeneous. 

To get a viable model we have to choose param- 
eters so the energy density in the relativistic fam- 
ily is small enough to avoid spoiling light element 
production. The typical energy CeS of a relativistic 
second family particle is dominated by its momen- 
tum, which scales with the expansion of the uni- 
verse as a{t)~^, because the expansion stretches 
the de Broglie wavelengths. Thus, as long as the 
second family remains relativistic, its mean energy 
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density is ps ~ ris^cS a{t)~'^, as usual for rel- 
ativistic particles. The ratio of the mean energy 
densities in the two DM families is 



(41) 



The first factor in the last expression must be 
larger than unity, say by a factor of 100, to keep 
(j) locked to the second family even as density con- 
centrations in the first family develop. To avoid 
affecting the standard model for the origin of the 
light elements we want the energy density in the 
rclativistic second family to be small compared to 
the thermal background radiation, meaning the 
present density is {ps/p)o ^ 10~*. Thus initial 
conditions for the particle momenta must be such 
that Ceff ~ kes/ao < I0~'^ys(l)*. 

Under the above assumptions, and neglecting 
dV/d(j) for the moment, the DE field equation 
when (h is near zero is 



V^ct)/a^ = t/>s|(/)|/eeff - yn{l + S). 



(42) 



In the first source term we have written the rcla- 
tivistic correction (eq. [14]) in terms of the effec- 
tive mean particle energy, e^s, as in equation (23): 



(Vl^) =2/s|0(x)|/6eff. 



(43) 



We have dropped the time derivatives of (j) because 
we are assuming the field is locked to a value near 
zero. The number density in the second (rcla- 
tivistic) family is nearly independent of position. 
The density contrast in the nonrelativistic DE is 
(5(x, t). On scales small compared to the Hub- 
ble length the second family particles see a nearly 
static potential, so the particle energy is con- 
served and thus independent of position along its 
trajectory. 

The space average of equation (42), which ne- 
glects dV/ d(j), gives the mean field value. 



<Pb = €csyn/{ylns). 



(44) 



This relation in equation (43) reproduces the con- 
dition that the mean inverse Lorentz factor for the 
second family satisfies 



(a/i - = yn/{ysns) < 1. 



(45) 



The departure from the mean of equation (42) 
(with (pi, > since we are taking dV/dcj) = for 



the moment) is 

V^(?ii/a^ = ylns(t>\/ees - ynS. (46) 
The Fourier transform is 

yn5{k) 



L(k) 



(47) 



The Green's function thus has a Yukawa form, oc 
exp ~r/r^, with cutoff length r., = y^i^Jff/yfn^. 
This scales with time as oc a{t), so the comoving 
cutoff length is constant. 

If the Hubble parameter H is dominated by the 
mass density j/n^* in the nonrelativistic DM, the 
cutoff length satisfies 



2 GecS(t>*yn I Ps ( yn 



P P \ysn 



■ (48) 



As argued above, the last two factors are at most 
10-* and (10-2)2, so that 



(^^^5)0 < 10-^0-'/^. 



(49) 



To produce r^iao) of order the Hubble radius 
would require f3 ~ 10-* and thus 0* ~ lO^rnp;, 
which is disagreeably large on theoretical grounds. 
On the other hand, (3=1 implies r5(ao) < 1 
Mpc which is well within the non-linear cluster- 
ing length, and so requires closer analysis to test. 

5.2. Late Time Transition 

The behavior of (f) when the number densities 
become small enough to allow the field to evolve 
depends on the self-interaction potential V{(j)). In 
linear perturbation theory for the field, the condi- 
tion that the second family particle number den- 
sity is large enough to hold </> constant and close 
to zero is 



dV/d(p = ±y,n,(v/l - v^) + yn{l + 6). (50) 

This assumes dV/d(j) > at ^ ~ 0. The negative 
sign in the first term on the right hand side applies 
when (/) > 0, and the positive sign when the par- 
ticle number densities are small enough to allow 
dV/d(j) to pull to a slightly negative value. Be- 
cause the second family is rclativistic, is nearly 
homogeneous, as we have discussed. This means 
the reciprocal Lorentz factor (eq. [43]) must be 
a function of position, balancing the irregular DM 
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mass distribution in the first family represented by 
the number density contrast (5(x, t). We will check 
below that even when the mass density contrast is 
non-linear, equation (50) is a good approximation 
for the perturbation to the field. 

We first consider the case where dV/dcj) can be 
neglected, so equation (50) is 



Vsn, 



=yn(l + 5). 



(51) 



There comes a time when the value (5max of the 
density contrast within the strongest concentra- 
tions of the massive DM family is large enough to 
satisfy j/sris = yn{l + 5max)- This forces the sec- 
ond family to become nonrelativistic in the neigh- 
borhood of (5max- Further expansion of the uni- 
verse increases the density contrasts, causing </> in 
the vicinity of a DM mass concentration to in- 
crease. When ys<l> > Ceff) second family particles 
are pushed out of the regions of first family con- 
centrations because their energy is not sufficient 
to allow them to have such a large mass. If this 
rearrangement is happening on length scales much 
smaller than the Hubble length, the DE field equa- 
tion is dominated by the spatial derivatives, and 
we have 



j/,n,(x,t)(v/T^) -yn(l + ,5), (52) 



for the case that dV/dcj) can be neglected. As in- 
dicated, the second family number density Ug is 
now a fimction of position, because less energetic 
particles are excluded from the concentrations of 
the first family. The spatial mean of equation (52) 
says that the inverse Lorentz factor averaged over 
all second family particles satisfies 



(53) 



Since we arc assuming the ratio on the right hand 
side is smaller than unity, the pools of relativistic 
second family particles in the regions between the 
concentrations of nonrelativistic DM are always 
able to hold the field value close to zero in the 
voids between the concentrations of galaxies. 

We have been assuming the field value within 
a concentration of the first family is less than 
so these particles are nonrelativistic. To estimate 
the extent of the mass shift of first family particles, 
consider a concentration of A'' ~ fiE? nonrelativis- 
tic DM particles drawn from an initially homoge- 
neous patch of size R into a concentration with 



density contrast 5 over a size r. The typical shift 
in the DE field value averaged over this concentra- 
tion is (^r ~ yN/r. Assuming the expansion rate 
is mainly due to the mass density in this family so 
the Hubble parameter satisfies H"^ ~ Gyncj)*:, we 
get 



(54) 



The fractional shift in the scalar field value is 
largest in the largest mass concentrations. For 
example, the density contrast in a large galaxy 
is about 6 = 10® at r = 10 kpc, which gives 
(pr/tp* ~ 10~^/3, and in a rich cluster at r = 2 
Mpc where S ~ 100, <i)r/<P* ~ lO"''/?. That is, 
as long as /? is not large the first family particle 
masses are only slightly perturbed by the spatial 
variation of (j), and the force on massive DM par- 
ticles from the gradient of </> is well approximated 
as /3 times the gravitational attraction (eq. [18]). 

Finally, let us briefly consider what happens 
when dV/ dip in equation (50) is positive and large 
enough to pull <p to negative values. This makes 
the second family first become nonrelativistic in 
the voids, where the first family number density 
is low. When this happens the field in the voids 
moves toward the minimum oi V , at (p = —(pg. If 
ys4>s ^ Ceflf then the DE field pushes the second 
family out of the voids. Since the mean particle 
number densities are decreasing as the universe 
expands, the potential V must eventually pull the 
DE field away from zero everywhere, producing a 
second family of nonrelativistic dark matter. 

When the DE field is no longer locked to the 
zero of mass for the second family there is a fifth 
force, but it need not parallel the peculiar gravi- 
tational attraction of the dark matter because the 
space distributions of the two families may differ. 

A model in which the lock on the DE field has 
broken before the present epoch and produced a 
large fifth force within concentrations of galaxies 
would not be acceptable, but a model with a large 
fifth force in the voids between the concentrations 
of large galaxies might be quite interesting, as a 
way to understand why the voids are so empty. 

6. Concluding Remarks 

We have explored physical processes in models 
where the DM particles have a Yukawa coupling 
to a scalar field that can be the source of the DE. 
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Our central conclusion is that parameters and ini- 
tial conditions in such models can be chosen so 
that the model is viable but significantly different 
from the standard ACDM cosmology. It might be 
useful to follow this up by considering whether the 
interesting range of initial field values, just some- 
what larger than the Planck mass for the models 
in §4, could naturally follow from models for the 
very early universe. It would also be useful to 
know whether the number density of DM parti- 
cles required to give the observed value of can 
be naturally understood, for example by gravita- 
tional production. 

As illustrated in §4 for a single DM family, the 
constraints on cosmological parameters derived by 
fitting the model for the dark sector to the ob- 
servations can differ from what is obtained from 
fitting to the standard model. Such alternative 
models are therefore useful as foils to ACDM, for 
the purpose of evaluating the empirical constraints 
on the cosmological parameters 

We have also considered two-family models, 
that can lead to a rich and interesting cosmology. 
As discussed in §5, initial conditions can be chosen 
so the DE field is locked to the zero of mass for the 
more numerous family, and remains so up to the 
present epoch. This removes the evolution of the 
DE field and the evolution of the mass of the non- 
relativistic DM particles, and it can suppress the 
fifth force in the dark sector. This is an example of 
how a dynamical model for the DE can be observa- 
tionally indistinguishable from Einstein's cosmo- 
logical constant. The situation changes when the 
DE particle number density becomes low enough 
to free the DE field. When this happens the be- 
havior can become quite different from the stan- 
dard model. 

Models of the type studied here, in which one 
or possibly several families of dark matter have 
masses set by their interaction with a dynamical 
scalar field, arc a useful cautionary example of an- 
other point; the empirical evidence on how the 
universe has been evolving up to now may be a 
rather deceptive guide to its physics or to its fu- 
ture. 
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